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COUPLING ALL THE LEVY STOCHASTIC AREAS OE 
MULTIDIMENSIONAL BROWNIAN MOTION 

By Wilfrid S. Kendall 
University of Warwick 

It is shown how to construct a successful co-adapted coupling of 
two copies of an n-dimensional Brownian motion {Bi,... ,Bn) while 
simultaneously coupling all corresponding copies of Levy stochas¬ 
tic areas f BidBj — f BjdBi. It is conjectured that successful co¬ 
adapted couplings still exist when the Levy stochastic areas are re¬ 
placed by a finite set of multiply iterated path- and time-integrals, 
subject to algebraic compatibility of the initial conditions. 


1. Introduction. A probabilistic coupling of two random processes is a 
construction of both processes on the same probability space, building in use¬ 
ful dependence between the two processes. This paper discusses couplings 
of two Markov processes with the same law of evolution, begun at different 
points, and constructed so as to join together (to couple) at some random 
time; the coupling is said to be successful if the two processes couple within 
finite time almost surely. There are other kinds of couplings relating to 
monotonicity, or to approximation; successful couplings are useful for prob¬ 
abilistic gradient estimates, for studying the rate of convergence to statistical 
equilibrium, for relating behavior of random processes to the geometry of 
the state-space, and (in more developed formulations) as a key component 
in perfect simulation algorithms. The present paper focuses on a particular 
question to do with coupling constructions for Euclidean Brownian motions: 
namely, whether one can couple successfully not only the Brownian motions 
themselves, but also sets of path functionals. We shall show that one can 
couple successfully not only two copies of a Brownian motion (Hi,..., Bn), 
but at the same time all the corresponding pairs of Levy stochastic areas 
/ Bi dBj — f Bj dBi of the two copies. This appears quite remarkable to the 
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author: one is able to couple so much despite controlling only the correlations 
between the two copies of the Brownian motion. 

Extensive treatments of probabilistic coupling can be found in [16] and 
[22], so a short summary of the relevant history will suffice. Lindvall [15] 
was the first to consider coupling for Brownian motion; he described the 
classic reflection coupling (couple two n-dimensional Brownian motions by 
making one of them to be the reflection of the other until they meet). This 
was followed up by Lindvall and Rogers [17], who discussed generalizations 
to the case of diffusions. There is a significant distinction to be drawn here. 
It is typically very much easier to find explicit descriptions of couplings 
when the two processes in question are required both to be co-adapted to 
the same hltration of u-algebras, in particular when the driving Brownian 
motions have increments which are independent of their common past. In 
[17] (and throughout the present paper) the search is for co-adapted cou¬ 
plings, and therefore stochastic calculus can be used to provide very explicit 
descriptions. 

Ben Arous, Cranston and Kendall [1] were the first to consider the possi¬ 
bility of what one might call exotic couplings, in which one seeks to couple 
co-adaptively and simultaneously certain path functionals as well as the pro¬ 
cesses. They described co-adapted couplings for the single stochastic area of 
planar Brownian motion, and also for the time-integral of scalar Brownian 
motion J B dt. Price [20] showed in her thesis how to extend the second case 
to couple the twice-iterated time-integral / / Bdsdt, and Kendall and Price 
[13] used a different method to show the existence of a successful co-adapted 
coupling for B and any finite set of iterated time-integrals f ■ ■ ■ J B ds ■ ■ ■ dt. 
The present paper continues this theme by extending the first result of [1] 
to n-dimensional Brownian motion and all possible Levy stochastic areas 
/ BidBj — BjdBi. It now seems reasonable to formulate a general conjec¬ 
ture that successful co-adapted exotic coupling is possible for any finite 
combination of multiply iterated path- and time-integrals of Brownian mo¬ 
tion (for compatible initial conditions), though it is clear that new ideas will 
be required to make further progress. The theory of Lie group symmetries 
of stochastic differential equations supports the expectation that resolution 
of the general conjecture would lead quickly to coupling constructions for 
wide classes of hypoelliptic diffusions. 

At present the main motivation for studying exotic coupling lies in the im¬ 
portance of coupling as a general concept, and the consequent desirability of 
understanding how far one can go in coupling large sets of path functionals. 
However, it does seem not unreasonable to hope for future useful interac¬ 
tions with rough path theory [18], where stochastic areas play a central role, 
and conceivably also for its use in lifting restrictions on the new methods of 
exact simulation of stochastic differential equations [2]. 
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It should be noted that there is significant theory concerning nonco- 
adapted couplings. If the co-adapted constraint is lifted, then one may 
construct maximal couplings [7, 8, 19] which couple at the maximum rate 
permitted by the total variation bound of the coupling inequality. These 
couplings have strong relationships with potential theory, and will in gen¬ 
eral be hard to construct (but see the striking results of Rogers [21] on 
random walks). Hairer [9] used a restricted kind of nonco-adapted coupling 
at time oo to study hypoelliptic diffusions, Hayes and Vigoda [10] used fi¬ 
nite look-ahead couplings to gain definite improvements on coupling rate 
in an application to randomized algorithms, while Burdzy and Kendall [4] 
studied the cost of the co-adapted property. In our case it is a simple mat¬ 
ter to demonstrate the possibility of successful nonco-adapted coupling of 
Levy stochastic areas as a consequence of Hormander regularity of the cor¬ 
responding n -|- n(n — 1)/2-dimensional hypoelliptic diffusion. The point of 
the present paper is to deliver an explicit successful co-adapted coupling 
construction; the existence of this is not implied by regularity theory. 

The paper is organized as follows. Section 2 addresses some general con¬ 
siderations related to stochastic control, which help to focus the problem 
on specihc coupling strategies and to introduce notation. Section 3 gives a 
new approach to the two-dimensional problem treated by [1]; this prepares 
the way for the main results of the paper which are stated and proved in 
Section 4: namely, that successful co-adapted coupling is achievable for n- 
dimensional Brownian motion and its n(n — l)/2 associated Levy stochastic 
areas. The concluding Section 5 considers a couple of complementary issues, 
and formulates a general question concerning coupling of sets of iterated 
path integrals, for which the answer is conjectured to be in the affirmative. 

2. Coupling, control and convexity. It is helpful to bear in mind a stochas¬ 
tic control-theoretic perspective for coupling problems concerning co-adapted 
Brownian motions (see Borkar [3] for a useful survey on stochastic con¬ 
trol; Chen [5] elicits some connections between control and coupling, while 
Jansons and Metcalfe [12] carry out some numerical investigations). As re¬ 
marked above, a co-adapted coupling of two n-dimensional Brownian mo¬ 
tions A and B means that A and B are both adapted as Brownian motions 
to the same hltration of cr-algebras t > 0} (thus in particular both incre¬ 
ments At-es ~ As and B^+s — B^ are independent of ^t)- The most general 
co-adapted coupling can be specified using Ito stochastic calculus: 

(1) dA = l^dB+^ dC, 

where J, J are predictable (n x n)-matrix-valued processes, and C is a fur¬ 
ther n-dimensional Brownian motion adapted as a Brownian motion to the 
filtration {^t : t > 0} and independent of B. 
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Thus the coupling is specified by giving a control in the form of a pair of 
predictable matrix-valued processes J, J. These must satisfy certain condi¬ 
tions if (1) is indeed to define an n-dimensional Brownian motion A: mul¬ 
tiplying stochastic differentials to obtain differentials of quadratic variation 
(following [11]), and bearing in mind the independence of B and C, it fol¬ 
lows that A is a Brownian motion if and only if the following matrix-valued 
random measure equation is satisfied: 

ldt = (dA) X {dAf 

(2) = (I^(dB) X (dB)^I) + (l’^(dC) x (dC)^I) 

= + f i)dT 

where I is the n x n identity matrix. 

The matrix process J expresses the infinitesimal correlation of C with 
B: from (2) it follows that such matrix processes are characterized by lying 
(almost) always in the convex compact set defined by 

(3) 0 < 1^1 < I, 

where 0 is the (ri x n) zero matrix, and the inequalities are to be interpreted 
using the usual spectrally based partial ordering for symmetric matrices. 
An application of the Cauchy-Schwarz inequality to v^J^J v shows that 
the set of extreme points of this compact convex set can be identified as the 
topological group of orthogonal matrices 0{n). 

Our coupling problem will be solved by designing a predictable process 
J such that A and B couple at some finite random time simultaneously 
with all their stochastic areas (/ Aj dAj — f Aj dAi coupling with f Bi dBj — 
J BjdBi, etc.). Suppose that it is possible to arrange this in terms of a 
stochastic control problem which is regular enough to possess an objective 
function leading to a bounded value function B(t, A,B) (where perhaps t is 
replaced by some other additive functional such as time spent in a specified 
region). Being a value function, B(t, A,B) is a supermartingale in general 
and is a martingale exactly, when the control J is optimal. If I^(t, A,B) 
is appropriately smooth, then Ito’s formula may be applied. This together 
with the Brownian nature of A and B shows 

dV{t,A, B) = Vodt + Vf dA + dB 

+ ^ tr(Yi]^) dt+ ^ tr( V.^ .,) dt + tr(JV^„j dt 

(for a fixed orthonormal basis vi,... ,Vn, and first- and second-order deriva¬ 
tives Bq, Yi, V 2 , Y,,, Y 22 with dependence on t, A and B suppressed). 

Thus processes J which are optimal controls for such a regular problem must 
maximize 

tr(^i2). 
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which is linear in J. It follows that smoothness of an appropriate value 
function implies that optimal control processes J must (almost surely, for 
almost all time) lie in the region of extreme points of the convex compact 
region of controls, and so must satisfy the orthogonality condition 

(4) 

in brief, almost surely J(t) £ 0(n) for almost all times t (and hence J = 0). 

The impact of these considerations for our coupling problem is entirely 
heuristic, since we do not have any particular objective function in mind 
other than desiring to show that it is possible to couple Brownian motions 
together with their stochastic areas. (Indeed we will not even check that 
our resulting coupling strategy is admissible, in the sense of being optimal 
for some objective function: it is not a priori at all clear whether a success¬ 
ful coupling exists and therefore optimality with respect to some arbitrary 
objective function is of less value than conceptual simplicity!) The above 
remarks encourage a search for simple couplings among those which use 
0(n)-valued processes J to construct A = / in terms of B, without 

any need of further randomness from C. Since 0(n) has two topological 
components, made up of SO{n) and the coset of rotated reflections, it also 
follows that we should expect to consider coupling strategies which involve 
discontinuous transitions between one control and another; and this is in¬ 
deed what may be observed for the successful coupling strategies described 
in Theorem 4 (for the planar case) and Corollary 7 (for the general case) 
below. 

3. The planar case. We first review the planar case (dimension n = 2), 
which permits a simpler treatment than the general re-dimensional case but 
introduces most of the key ideas. The planar case was first dealt with in [1], 
using controls J not all lying in 0(2), though it was noted in passing that 
there was a possibility of coupling using only reflection and synchronous 
coupling (as defined in Definitions 2 and 3, J is a reflection matrix or is 
an identity matrix). Ben Arous and Lyons have shown in unpublished work 
how to implement reflection/synchronous coupling for the planar case in 
a rather direct way, which resembles the low-dimensional case of the Ben 
Arous et al./Kendall and Price treatment of Kolmogorov diffusions (time- 
integrals and twice-iterated time-integrals together with scalar Brownian 
motion). Here we show how reflection/synchronous coupling may be set up 
using simple and largely state-dependent coupling rules. 

First recall from [1] that it is sufficient to couple {Bi,B 2 ) and (^ 1 ,^ 2 ) 
together with the invariant difference of their stochastic areas, 

(5) 21 = j {AidA2- A2dAi) - j {BidB2- B2dBi) + A 1 B 2 - A 2 B 1 . 
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In fact 21 then has a geometric interpretation: it measures the stochastic 
area swept out by moving first along the A path, then linearly from the end 
of the A path to the end of the B path, and then back along the B path 
to its starting point. It turns out to be natural to think of 21 as the (1,2) 
coordinate of an antisymmetric matrix 


21 = 




Consider the summary quantities 


( 6 ) 


V = ^{Ai-BiY + {A2-B2)\ 
U = sgn(2t)^tr(^^^) = \/22l. 


These are semimartingales at least until one of them vanishes. Stochastic 
calculus can therefore be used to compute the stochastic differential drifts 
Drift dU and Drift dV (the differentials of the locally bounded variation com¬ 
ponents of the Doob~Meyer semimartingale decompositions of C/, V) and 
the products of differentials {dU)'^, {dV)"^ and {dV) x {dU) (the differen¬ 
tials of the corresponding quadratic variation and covariation processes). In 
doing this, it is convenient to define the quantities 5ii, S 22 and A 12 from 
the symmetrization and the antisymmetrization of the control J: working in 
orthonormal coordinates based on the vector A — B and its perpendicular. 


( 7 ) 




Ai2 \ 

0 )■ 


The results of these computations are summarized in the following lemma: 


Lemma 1. 

{dVf = 2{l- Sii)dt, Drift dC = 

(8) {dV) X {dU) = -2V2Ai2Vdt, 

{dUf = 4(1 + S 22 )V^ dt, Drift dU = 2 V 2 A 12 dt. 

Details of the calculations are left as an exercise for the reader, who may 
alternatively view them as a special case of the multidimensional case treated 
in detail in Lemma 5. 

Here are two important coupling strategies, defined by specifying the cor¬ 
responding control J. 

Definition 2. Reflection coupling is defined by choosing J to be the 
orthogonal matrix giving reflection in the line normal to the vector A — B: 
thus 5ii = — 1 , S 22 = I 5 Ai 2 = 0 in our chosen coordinate system. 
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Using Lemma 1, reflection coupling yields 

{dVf = 4.dt, Drift dU = 0 


(9) 


{dV) X {dU) = 0, 

{dUf = 8V^dt, Drift dC/ = 0 


so that V moves as a scalar Brownian motion at least until it hits 0, and U 
moves as a scalar Brownian motion subject to a U-dependent time-change. 

Definition 3. Synchronous coupling is defined by choosing J to be the 
identity matrix; thus 5ii = S 22 = 1, ^12 = 0. 

Using Lemma 1, synchronous coupling yields 


{dVf = 0, Drift dU = 0, 


( 10 ) 


{dV) X (dU) = 0, 

{dUf = dt, Drift dU = 0, 


so that V is held constant, while U continues to move as a scalar Brownian 
motion with rate dependent on V in the same way as for reflection coupling. 

Under both these strategies U and V remain semimartingales for all time. 

It is possible to derive these results for both couplings without making 
explicit use of stochastic calculus, simply by considering the geometry of the 
planar Brownian paths and their invariant difference of areas. 

The considerations of Section 2 suggest that if coupling is at all possible 
for the planar case using only symmetric J, then it should be achievable by 
combining these two controls, since (8) shows that the other two extreme 
controls (5ii = ±1, S 22 = ~1) lead to positive drifts for V without apparent 
gains for U. 

Since U scales as V^, and since it is desirable for coupling purposes to 
reduce the size of U if ever it gets large relative to V, it is natural to consider 
coupling strategies described loosely as follows: for fixed k > 0, 

while use reflection coupling; 

while use synchronous coupling. 

This involves a discontinuous change of regime as {U, V) crosses over the 
boundary = k^V^. The discussion in Section 2 has prepared us to expect 
such discontinuities. A precise description of a successful strategy of this 
kind is formulated in the following theorem, which is the principal result of 
this section. 

Theorem 4. Suppose that initially Uq = 0 but Vq > 0 (this can always 
be arranged by first using reflection coupling to make V positive, and then 
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using a session of synchronous coupling to reduce U to zero). Fix a small s > 
0; consider the control which alternates between reflection and synchronous 
couplings using ^''down-crossings'’'’: 

- ifU’^/'V'^ has not yet visited , then use reflection coupling; 

- if has attained the level {k — efl since most recently visiting n?, 

then use reflection coupling; 

- otherwise use synchronous coupling. 

This coupling is almost surely successful in finite time: {U,'V) visits (0,0) in 
finite time. 

Clearly one could consider the limiting case e ^ 0 and use local time and 
excursion theory; however, it turns out to be simpler to analyze the process 
as given. 

Proof of Theorem 4. Dehne the indicator process by = 1 
when either has not yet visited or has attained the level 

(k — efl since most recently visiting and otherwise set = 0. Then 
the coupling strategy prescribed in the theorem statement corresponds to 
the stochastic differential system 

{dVf = dt, Drift dV = 0, 

(11) {dV)x{dU)=0, 

{dUf = 8P2 dt, Drift dU = 0. 

This is solvable up to the time when U and V both vanish: one may piece 
together solutions of the smooth systems defined by (9) and (10). Under this 
stochastic differential system the process V evolves as a Brownian motion 
of rate 4 interrupted only when /'V^ makes down-crossings from nf to 
(k — e)^, and during these interruptions V is frozen. These down-crossings 
each take a finite amount of time, and only finitely many occur in bounded 
closed time intervals before U and V both vanish; consequently V either 
hits 0 at a finite time or converges to 0 at time oo. Since V is constant when 
jV^ > nf, continuity considerations show that U/V‘^ —> 0 as U —> 0, and 
therefore coupling must occur when V hits 0. Thus the crux of the matter 
is, will U —> 0 at a finite time? 

To analyze this question, apply Lamperti’s [14] observation (as used to 
great effect in [23], e.g.) to the stochastic differential system (11). Consider 
a random time-change under which K = log(U) behaves as an (interrupted) 
Brownian motion with constant negative drift. The time-change T{t) is de¬ 
fined by 

( 12 ) 


Adt = V^ dr. 
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Writing W = U/V'^, the stochastic system for K and W then follows by Ito’s 
formula: 

{dKf = dr, Drift dK = dr, 

(13) {dK) X (dW) = 2N^^^WdT, 

{dW)‘^ = 2{l + 2N^^^W‘^)dT, Drift diy = 3iV(^)W dr. 

It is required to show that elapsed t-time until K —oo (equivalently D = 0) 
is finite, which is equivalent to showing 

POO 

(14) / e^^dT<oo. 

Jo 

Since V diffuses as Brownian motion of rate 4 when = 1 and is otherwise 
frozen, it follows that the integral dr is a Brownian first-passage 

time and therefore is finite. Accordingly, it is enough to show 

POO 

(15) / (1 — dr < oo. 

Jo 

Let < cr^ be the start and finish times (in r-time-scale) of the nth 
down-crossing of = U'^/V^ from k? to (k — e)^. But Nr^'^ = 0 exactly 
when r lies in the union of the stopping-time intervals [(r®,(T^], so therefore 

POO ^ 

(16) / (l-lv(^))e^^dr = ^e^^<(o-^-0, 

•^0 n=l 


since V = remains constant for r G [(T^,cr^]. 

Conditional on Kfj^: n = 1,2,..., the durations are independent 

Brownian first-passage times of different rates. Consequently 


E 


(17) 


\ n=l 


K„3 : n = 1,2,.. 

^ T3. ' ' 




= exp — ^ X £ I, 
V n=l 


using the formula for the moment-generating function of a Brownian first- 
passage time. 

Consider now the times ul — (t{, — a^, - ■ ■ between successive down- 

crossings. These are independent, identically distributed, and of finite mean, 
since their common distribution is the r-time for the regular real-line diffu¬ 
sion W (with = 1) to hit one of ±k when started at k — e. Thus by the 
strong law of large numbers it follows that almost surely 

1 ” 

^ E[u| - Cj(] > 0 

fi 

m=l 
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(defining cJq = 0 ). 

But equally iiT is a Brownian motion with constant drift of — ^ on the 
interrupted r-time-scale / dr, and therefore almost surely 

_ ^ _ 

It follows that almost surely 

(18) i|l^_iE|,r|-<Tf]<0, 

and hence is almost surely finite. 

Consequently (17) shows that 


E 


exp 



2K„ 


n = l,2,. 


is almost surely positive, and so 


OO 



n=l 


has a positive chance of being finite, even when conditioned on dLo-s: n = 

1,2,_But the e '^"(<7^ — cr^) are independent under this conditioning, 

and so by the Kolmogorov zero-one law and (16) it follows that 

roo 

(19) / {1 - dr < OO 

JO 

with probability 1. It follows that coupling under this strategy almost surely 
succeeds after a finite time. □ 


Further development of this line of reasoning delivers an explicit construc¬ 
tion of the limiting case e —> 0 using local time and excursion theory, a single 
elliptic partial differential equation for the moment-generating function 

E[exp(-pr)] 

of the coupling time T for all p using scaling, and estimates for exceedance 
probabilities of the coupling time. We do not consider these topics here, but 
instead proceed to the multidimensional case. 

4. The n-dimensional case. The hrst step is to establish the stochastic 
differential system (6) for Euclidean separation and invariant difference of 
stochastic areas, working in general n-space (n > 2). First introduce new 
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coordinates based on X = A ~ B and Y = A + B, where A and B are co¬ 
adapted n-dimensional Brownian motions satisfying (1). Using ltd calcnlus 
for the vectors dX. and dY, 

dX dX^ = 2(1-g)dt, Drift dX = 0, 

(20) dYdX^ = 2Adt, 

dY dYj = 2{l + g)dt, Drift dY = 0, 

where S = ^(JT J^) and A = ^(J — J^) are the symmetrized and antisym¬ 
metrized matrices corresponding to J. 

Applying the ltd formula to = X^X (the square of the length of X), 
it follows that while V remains positive 

(21) (dV)'^ = 2{l-u^gE)dt, Drift dU = ~ ^ ~ ~ " =-^ dt, 

where u = X/V is a normalized configuration vector defined by X = A — B. 

Now consider the antisymmetric matrix % determined by invariant differ¬ 
ences of stochastic areas of the form of (5); 

21 = J{^i dAj — Aj dAi) — J{Bi dBj — Bj dBi) -|- AiBj — AjBi. 

Since A = ^ (Y -|- X) and B = ^ (Y — X), calculation shows 

(22) d'Qiij = Xi dYj — Xj dYi — 2Aij dt. 

Hence 


d2ljj X dOirs 

(23) = XiXr dYj dYs - XjXr dYi dYs - XiX^ dYj dYr -L XjX^ dYi dYr 
= 2{X,Xr a + ^),s - XjXr a + ^)« 

- XiXsil + g)jr+XjXsil + g)ir)dt. 

Setting U = tr(^^^), since 

(24) d{U^) = 2UdU + {dUf = d^ij + {d2iijf) 

* i 

it follows 

4C/2 {dUf = 4 X! Z! Z! Z! d^rs 

i j r s 

= 32 Z Z Z Z ^ij^rsXiXr (I + g)js dt 

i j r s 

= 32X^g^{I + §)gx dt = 32^3^ (I + dt. 


(25) 
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Here 3 =2i;/H is a normalized configuration matrix [with tr(3^3) = 1, anti¬ 
symmetric so 3^ = —3 and = 0]. The second line of (25) follows from 
the first by applying (23) and then exploiting the symmetry of I -|- S and 
the antisymmetry of 21. 

On the other hand, from (24), 

2 [/Drift dH 

= Drift d(D2) - {dU f 
= Drift^ ^(22lij d%ij + {d^ijf) - {dUf 

* 3 

= 4tr(2t^A) dt — {dU)'^ 

( 26 ) + ^^3 

* 3 

= AiT{^^ dt - {duf 

+ 2 E E(2^*' (I + Dii - 2^*^J (I + gh ) dt 

i 3 

= 4tr(§^A) dt - {dUf + 4(tr(I + g)- + g)]/)V^ dt 

= 4tr(^^A) dt -L 4(n — 1 -|- tr ^ — 2^^3^(| -|- ^)3l:)H^ dt, 

where the last line is derived from (25), evaluating trl = n, u^lu = 1. 
From(25) and (26) taken together, 

{dUf = dt, 

(27) Drift dU = 2 tr(3^A) dt 

T/2 

-L 2 (n — 1 -L tr^ — — 2 EE(| + ^)3k) -jj- dt. 

Finally, using the antisymmetry of 21, 

d{U^) d{V^) = WU dV dU 

=4j2x,dx.Y.i: 2tr-s d^rs 

i r s 

= -leEEAiiDEdt 

and so finally 

(28) dU dV = -4EEAii^ dt. 

Following the procedure of the planar case, now consider the behavior of 
K = log(F). As in Section 3, define W = U/V'^; however, we will consider 
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the behavior of K together with that of H = log(C/) rather than that of 
W = exp(i? — 2K). The next lemma follows from the calcnlations in this 
section so far. 


Lemma 5. For a general control J (with symmetric and antisymmetric 
components S and A), and defining a new (T-)time-scale by Adt = V‘^ dr as 
in Section 3, 

{dK)‘^ = dr, 

Drift dK = |(n — tr S — 2(1 — vJ'Sv)) dr, 

{dK) X {dH) = 

(29) {dHf = dr, 

Drift diL = ^ tr(3^A)--^ dr 

+ ^(re-l + tr^ - - 4i/'^3^(I + §)3l) dr. 

Proof. Use (21), (27) and (28), and Ito’s formula. □ 


The special cases of reflection and synchronous coupling now follow di¬ 
rectly. Reflection coupling is defined by 

(30) 

which implies 

^reflection ^ 4 = 0, 

trS = n —2, = —1. S3j£ = 3T; 

and consequently 

{dK)"^ = dr, Drift dK = —^dr, 

(31) {dK)x{dH)=0, 


dr 


dr 


(dH) =4 \\^m.\\ DriftdL7 = (n-l-411^1^11 ) 

Synchronous coupling is defined by 

^22^ jsynchronous _j 

which implies 

(g _jsynchronous _ q 

trS=ri, = 1. S3 ii = 3t, 


14 


W. S. KENDALL 


and consequently 

{dKf = id, Drift (iK = 0, 

(33) {dK) X {dH) = 0, 

{dHf = AW^vf DriftdF = (n - 1 - 4||^^f) 

Note that ||3l£|P is bounded above by 1/2, since the nonzero eigenvalues of 
the antisymmetric matrix 3 ah have multiplicity 2 and the sum of squared 

eigenvalues is tr(3^3) = 1- So if n > 3, then H is a nonconstant submartin¬ 
gale under both controls; it follows that there is no hope of coupling higher¬ 
dimensional stochastic areas by using only synchronous and reflection cou¬ 
pling. Instead we analyze the more complicated case of general orthogonal- 
matrix controls. 

Consider the case of a rotation coupling defined adaptively by a matrix 
exponential 

(34) l™*'^'“‘^(6»2)=exp(02). 

Here 3 is an antisymmetric matrix satisfying tr(3^3) = 1, so that J''°*®'*'°“(03) 
is indeed a rotation matrix, and moreover a finite Taylor series expansion 
produces an approximation which can be bounded: 

^ = cosh(03) = i(l™t"t““(02) + J™'"*““(-02)) = I - + 0^0(1), 

A = sinh(6»2) = i(A°*'^*““^(02) -1™*'^'““(-02)) =9g + 0^2(1). 

Here the 0(1) terms in the errors signify matrices which vary from line to 
line but which can be bounded uniformly in 9 and 3- Hence 

tr S = n — ^ 0^0(1), = 1 — ^||3t|P + 9^0{1), 

tr(^^A) = 9tv{^^g) + 9^0{l), 1/^^^At = 9qu,gE) + 9^0{l), 

T^3^(i + I)3T=2||lidf+0^0(1), 

where again the 0(1) terms in the errors (both here and in the following 
exposition) vary from line to line but are bounded uniformly in 0, 3 and 
the configuration matrix 3- For the sake of simplicity we choose 9 = —'ylW, 
3 = 3, and consider the ®ects of applying the adaptive rotational control 

fjT 'y‘^ 
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dr 

Drift dir = ^ + ^O(l) dr, 

(35) (dir) X (dii) = -iWlvf ^ + ^^0{l) dr, 

Driftdif = -(2 - {„ - 1 - 4||j£|p)) ^ + f^O(l)dT. 

The antisymmetric component of the control contributes a crucial 
term to the drift of H. This can be used to make H a supermartingale. 
(Incidentally, the choice ^ = ^ maximizes this particular term.) 

This motivates a direct construction of a successful coupling strategy, 
using a mixture of and 73 /kD) with adaptive choices of 

parameters. This delivers a positive chance of successful coupling for large 
initial values VTo of W: 

Theorem 6. Consider the adaptively mixed coupling 

T ^ jreflection \ f -\ ^ ^ jrotation f 

= “IT 2 = +'v w^)= 1 

(36) d = = 2(^^lK + ^(1 - 2|||i2||2)) , 

7 = 7(^,t) = 2{hh + n - 1 - 4||^^|p), 
defined so long as 

> do = + ihH + n— 1)^. 

This coupling strategy has a positive probability of being successful within 
finite time if Wq > , where Wq is the initial value of W at time 0 and 

is a certain finite threshold defined by (39) below, so long as we choose 

(37) t)< UK < Th < hK- 

Moreover, the coupling strategy will succeed almost surely if W stays above 
the threshold for all time. 

Recall from the discussion after (33) that ||3 t|P is bounded above by 1/2. 
So d — 2pK as given above is always nonnegative (as is 7 — 2fiH if n > 3). 


Proof of Theorem 6 . The effect of the mixed control can be evalu¬ 
ated as a convex combination of the systems of reflection coupling (30) and 
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rotation coupling (35): 


(38) 


rlr 

{dKf = (2/zx + {m+n-l- 4||1i£|P)2(1 - W^uf)) ^ 

Drift dK = -MK:^ + ^ dr, 

{dK) X {dH) = 2{nH + n-l-A\\^uf)\\^uf ^ + ^^dT, 


{dHf = 4||3ld 


dr 


0 (1) ^ 

+ ^h7-dT, 


= " 

Drift dH = + ^ dT. 

The 0(1) terms here may be taken to be bounded uniformly in the configu¬ 
ration vector u and matrix 3, and in W. Choose e so that 2fiK — > e > 0 

and set df = dr/W'^, and use the bounds on the 0(1) terms to define 
< oo as the smallest level w such that 


(39) 


Drift dK 


dr 


+ IJ^K 


< 


3’ 


Drift did 


dr 


+ I^H 


< 


3’ 


whenever IT > u;. Recall that In(lT) = H — 2K, so the calculations of (38) 
show that (dlnVT )^/dr is bounded, while 


(40) 


Drift d In IT 


dr 


- (2/Ux - hh ) 


<e 


whenever IT > . Now e was chosen so that 2jiK — > e > 0, so it follows 

by consideration of the law of the iterated logarithm that if initially VTo > 
then there is a positive chance that IT > for all time; moreover, 
this probability increases to 1 as ITo increases. In case IT > for all 
time, IT will grow at least linearly with rate — e > 0, and hence 

(by considering for progressively smaller e) 


(41) 


In IT 


2flK - 


as r —> oo. 

On this event of linear growth of IT > the approximations in (38) 
improve with time. Thus as r —> oo, so the same application of the law of 
the iterated logarithm leads to 


K 


H 


(42) 




-IJ'H- 
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In summary, there is a positive probability of both (41) and (42) holding 
so long as Wq > is sufficiently large; indeed this probability increases 
to 1 as Wq —> oo. If fiK and hh are both positive, then this ensures that 
V = exp{H) and U = exp(iC) both hit zero (delivering coupling of both 
position and all stochastic areas) at r = oo. 

In principle the coupling might still happen at t-time oo, in which case it 
would not succeed at finite time. However, 

(43) dT = ^=(^'^ dt = exp{2{K - H)) dt 
and therefore the coupling will occur at t-time 

poo 

(44) / exp{—2{K — H))df. 

Jo 

This will be hnite on the event of linear growth of W if the positive px and 
Hh are chosen not only to ensure that (42) holds but also so that px < fJ-H- 

Consequently there is positive probability of coupling occurring at finite 
time so long as we have arranged for fix and fin to satisfy (37). □ 

Corollary 7. The adaptive mixed coupling of Theorem 6 can be mod¬ 
ified by adding a synchronous coupling regime so as to ensure successful 
coupling in finite time with probability 1. 

Proof. We have already dealt with n = 2 in Section 3 above. If n > 3 
and if W falls below , so that the above procedure breaks down, then we 
can revert to pure synchronous coupling (33) [which holds K constant and 
allows H to evolve as a nonconstant submartingale as noted after (33)] until 
W does exceed , and restart the procedure. Consequently the above can 
be converted into a strategy which produces coupling at finite time almost 
surely. □ 

The coupling strategy described in Corollary 7 involves discontinuous 
transitions between synchronous and mixture strategies, fulfilling the ex¬ 
pectations of the heuristics at the end of Section 2. Provided we resort to 
time-dependent strategies, we can of course replace the mixed strategy by a 
time-dependent variation between reflection and rotation strategies; hence 
coupling can be achieved using only orthogonal controls. 

5. Complements and conclnsion. It is natural to ask whether anything 
might be gained by considering the full coset of coupling strategies alternate 
to the rotation strategies: what we might call the rotated reflection couplings 

( 45 ) ^rot-refl^ _ 2 ^T^) exp(02). 
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Applying the same reasoning as led to (35), we find that J™* 
the following effect: 





Dr 




This analysis would lead to a rather transparent coupling strategy if we 
could ensure that H was always a supermartingale under a suitable rotated 
reflection coupling for small "f/W; however, this is not possible for n> 3 
since it can be shown that 



(47) 


for = (I~TT^)3(I —TT^) with the maximum being achieved when 3 = 3^- 
This maximum vanishes when 3 is of rank 2 and i^\s a. nonzero eigenvector 
of 3, so the evolution of the configuration (^, 3) unavoidably affects whether 
or not the drift of H is negative. 

It is also natural to ask whether a more direct analysis can be made 
using the Carnot-Caratheodory distance for the relevant nilpotent Lie group. 
Recall that the Carnot-Caratheodory distance between the origin 0 and a 
point X with specified stochastic areas 21 is obtained by minimizing the 
Euclidean length of paths from 0 to x which produce the specified matrix of 
stochastic areas. A variational analysis shows that in general these paths are 
Cartesian products of circular arcs. A direct but laborious computation can 
be made of the stochastic calculus for the Carnot-Caratheodory distance in 
the two-dimensional case; unfortunately no useful picture seems to emerge 
from these computations. 

There are various further questions to be addressed about stochastic area 
couplings. Certainly it is possible to use the methods described here to derive 
estimates on coupling rates; these are not pursued for reasons of space and 
also because there is a much more substantial open question: 
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Can one co-adaptively couple not just the Brownian motions and their stochas¬ 
tic areas, but also all possible iterated path and time-integrals up to a fixed 
order of iteration? 

Here of course it is necessary to suppose compatibility of the initial con¬ 
ditions, to avoid obstructions caused by algebraic relationships between 
the various iterated integrals (see, e.g., the algebraic remarks of Gaines 
[6]). Kendall and Price [13] answer this question affirmatively for the one¬ 
dimensional case by using an implicit approach; the work of this paper shows 
that all singly iterated path-integrals can be coupled co-adaptively, since 
these can all be expressed as linear combinations of Levy stochastic areas 
and quadratic functions of Brownian coordinates. The general n-dimensional 
case is much more involved. We conjecture nevertheless that there is an affir¬ 
mative answer to the full multidimensional question given above. However, 
it is clear that new approaches will have to be tried here as in the one¬ 
dimensional case: the structure which facilitates the matrix-based approach 
of Section 4 is no longer available for higher-order iterated integrals. 
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